Majorana modes at the vortex ends of superconducting doped topological insulators 
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Recent experiments have observed bulk superconductivity in doped topological insulators. Here 
we ask whether vortex Majorana zero modes, previously predicted to occur when superconductivity 
is induced on the surface of topological insulators, survive even in these doped systems with metallic 
normal states. Assuming inversion symmetry, we find that Majorana zero modes indeed appear but 
only below a critical doping. The critical doping is associated with a topological phase transition 
of the vortex line, where it supports gapless excitations along its length. The critical point depends 
only on the orientation of the vortex line, and a Berry phase property, the SU(2) Berry phase 
of the Fermi surface in the metallic normal state. By calculating this phase for available band 
structures we determine that materials candidates like p-doped Bi 2 Te3 under pressure supports 
vortex end Majorana modes. Surprisingly, even superconductors derived from topologically trivial 
band structures can support Majorana modes, providing a promising route to realizing them. 



Majorana fermions, defined as fermions that are their 
own anti-particles unlike conventional Dirac fermions 
such as electrons, have long been sought by high energy 
physicists, but so far in vain. Of late, the search for Majo- 
rana fermions has remarkably shifted to condensed mat- 
ter systems [1—3], especially, to superconductors (SCs), 
where states appear in conjugate pairs with equal and 
opposite energies. Then, a single state at zero energy is 
its own conjugate and hence, a Majorana state or a Ma- 
jorana zero mode (MZM). These states are immune to 
local noise and hence, considered strong candidates for 
storing quantum information and performing fault tol- 
erant quantum computation [4]. Moreover, they show 
non-Abelian rather than Bose or Fermi statistics which 
leads to a number of extraordinary phenomenon [5] . 

Despite many proposals direct experimental evidence 
for a MZM is still lacking. While initial proposals in- 
volved the v = 5/2 quantum Hall state and SCs with 
unconventional pairing such as p x + ip y [3], a recent 
breakthrough occurred with the discovery of topologi- 
cal insulators (TIs) [6], which feature metallic surface 
bands. When a conventional s-wave SC is brought near 
this metallic surface, a single MZM is trapped in the 
vortex core [7]. Since then, several TIs were found to 
exhibit bulk superconductivity on doping [8] or under 
pressure [9, 10]. The normal phase of these SCs is now 
metallic, which raises the question: can a SC vortex host 
a surface MZM even when the bulk is not insulating? 

In this Letter, we answer this question in the affir- 
mative and in the process, discover a convenient way to 
obtain a MZM, which allows us to conclude that some ex- 
isting experimental systems should possess these states. 
Our proposal involves simply passing a magnetic field 
through a Tl-based SC, such as superconducting Bi 2 Te 3 , 
in which the doping is below a certain threshold value. 
We also find general criteria for SCs to host vortex MZMs 
and show that some non-TI-based SCs satisfy them too. 

A heuristic rule often applied to answer the above ques- 
tion is to examine whether the normal state bulk Fermi 



surface (FS) is well separated from surface states in the 
Brillouin zone. If it is, MZMs are assumed to persist in 
the bulk SC. While this may indicate the presence of low 
energy states, it is not a topological criterion since it de- 
pends on non-universal details of surface band structure, 
and cannot signal the presence of true MZMs. For MZMs 
to disappear, a gapless channel must open that allows 
pairs to approach each other and annihilate. We there- 
fore search for and offer a bulk rather than surface cri- 
terion. In this process, we have uncovered the following 
interesting facts. We assume inversion (I) and time re- 
versal (T) symmetric band structures, and weak pairing, 
since these lead to a technical simplification and capture 
many real systems, (i) The appearance of surface MZMs 
is tied to the topological state of the vortex, viewed as 
a ID topological SC. The critical point at which they 
disappear is linked to a vortex phase transition (VPT) 
where this topology changes. If verified, this may be the 
first instance of a phase transition inside a topological 
defect, (ii) The topological state of a vortex depends in 
general on its orientation. (Hi) Symmetry dictates that 
the normal state FS is doubly degenerate, leading to an 
SU(2) non-Abelian Berry phase [11] for closed curves, 
which determines the condition for quantum criticality of 
the vortex. This is a rare example of a non-Abelian Berry 
phase directly influencing measurable physical properties 
of an electronic system. Spin-orbit coupling is essential 
to obtaining the SU(2) Berry Phase, (iv) Using this cri- 
terion and available band structures we find that MZMs 
occur in p-doped superconducting Bi 2 Te 3 [10] and in Cu- 
doped Bi 2 Se 3 [8] if the vortex is sufficiently tilted off the 
c-axis. C-axis vortices in Cu-doped Bi 2 Se 3 are predicted 
to be near the topological transition. 

The Problem: Consider a 3D insulating band structure 
H, which we dope by changing the chemical potential /i 
away from the middle of the band gap. Now, add conven- 
tional 's-wave' even parity pairing A (in contrast to the 
p-wave pairing of Ref. [12]) and introduce a single vortex 
line into the pairing function A(r), stretching between 
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the top and bottom surfaces. We neglect the effect of 
the magnetic field used to generate the vortices, assum- 
ing extreme type II limit. When H is a strong TI, and /i 
is in the band gap, the pair potential primarily induces 
superconductivity on the surface states. In this limit it 
is known [7] that MZMs appear on the surface, in the 
vortex core. [13] Now consider tuning /i deep into the 
bulk bands. By modifying states well below fj,, one could 
tune the band structure to one with uninverted bands. 
One now expects 'normal' behavior, and the absence of 
MZMs. Therefore, a quantum phase transition must oc- 
cur between these limits and [i — \i c . 

To understand the nature of the transition, we recall 
some basic facts of vortex electronic structure, which 
are also derived below. Once \i enters the bulk bands, 
low energy Caroli-de Gennes-Matricon excitations ap- 
pear, bound to the vortex line [14]. These excitations are 
still typically gapped, although by a small energy scale, 
the 'mini-gap': S ~ A/(fcp£), (where kp is the Fermi 
wave- vector and £ is the coherence length. In the weak 
pairing limit kp£ >• 1). This small energy scale arises 
because the gap vanishes in the vortex core leading to a 
droplet of normal fluid, which is eventually gapped by the 
finite vortex size. However, the presence of the minigap 
is important, since it blocks the tunneling of the surface 
MZMs into the vortex line, and confines them near the 
surface. The closing of the mini-gap allows the surface 
MZMs to tunnel along the vortex line annihilate. 

Vortex as a ID topological SC: The VPT may be 
viewed as a change in the topology of the electronic struc- 
ture of the vortex line. The relevant energy scale is of 
the order of the mini-gap S <C A, with excitations lo- 
calized within the ID vortex core. The vortex admits 
particle-hole symmetry (C) but breaks T-symmetry and 
hence, belongs to class D of the Altland-Zirnbauer clas- 
sification [15]. Thus, the problem reduces to classifying 
gapped phases in ID within the symmetry class D, which 
are known to be distinguished by a Z 2 topological invari- 
ant [1]. The two kinds of phases differ in whether they 
support MZMs at their ends. The topologically nontriv- 
ial phase does and hence, corresponds to the fi < fi c 
phase of the vortex line. On raising /i, the vortex line 
transitions into the trivial phase, via a quantum critical 
point at which it is gapless along its length. This is remi- 
niscent of recent proposals to generate MZMs at the ends 
of superconducting quantum wires [16]. Note, since there 
is no 'local' gap in the vortex core, the powerful defect 
topology classification of [17] cannot be applied. 

The Hamiltonian is H = ^ Efc *fe^fc dG *fc where 
= (<4 t , cj,^, c_/4, -c-kt) and cJ, CT is assumed to 
have a — 1 ... N orbital components ct and 



where scalars like /i and A multiply the identity matrix 



l2JVx2JV- The band Hamiltonian, Hf,, is a 2N x 27V ma- 
trix with T symmetry: <3ytl*_ k o-y = H^, where a y acts on 
the spin, which yields the Hamiltonian structure above. 
When I-symmetry is also present, Hf- will be doubly 
degenerate, since the combined operation TI leads to 
a Kramers pair at every momentum. "H]? dG has parti- 
cle hole symmetry implemented by the transformation 
C = TlyfjyJC, where n matrices act on Nambu particle- 
hole indices, and JC is complex conjugation. A vortex 
given by A(r) = |A(r)|e -je , breaks T but preserves C. 

Role of vortex orientation: Consider a straight vortex 
along z. The dispersion £{k z ) in general, has a minigap 
5 as in Fig. 1. A topological phase transition requires 
closing of the minigap and reopening with inverted sign. 
Since the Z 2 topological index is only changed by an odd 
number of such band crossings, the only relevant mo- 
menta are k z = 0, it. Band touchings at other k z points 
occur in pairs at ±k z which do not alter the Z 2 index [1]. 
In the weak pairing limit, one expects the critical point 
fj, c to be determined by a FS property, which will be out- 
lined in detail below. Here we simply observe that the 
relevant FSs to consider lie in the k z — 0, n planes, the 
planes determined by the vortex orientation. This im- 
plies that the topological phase of the vortex, and hence 
\x c depend in general on its orientation. 

VPT in a lattice model: Before discussing the gen- 
eral criterion for a VPT, we present numerical and ana- 
lytical evidence in a specific lattice model from Ref. [18]. 
While the numerics explicitly demonstrate the phase 
transition, the analytical treatment of the continuum 
limit allows us to conjecture a Berry phase condition 
for the transition, which is later proved. The model 
is on a simple cubic lattice with two orbitals per site: 
ifjf u = r x dk • cr + TOfcT z — n where r, (<7j) are Pauli ma- 
trices in the orbital (spin) basis, d l k = 2t sin fcj, m,k = 

(M + mo J2i cos &*)) * = x i V-> z ; an d t, m o and M are pa- 
rameters of the model and (i is the chemical potential. 
The model is in the strong TI phase if— 3 < ^ < — 1. 
We add a mean field s-wave pairing to this Hamiltonian, 
insert a unit winding into the pairing function and diago- 
nalize the Hamiltonian numerically. We focus on k z = 0. 

Numerical results: Fig. 1 illustrates the evolution of 
the bulk vortex bound states, the dispersion within the 
vortex and the surface MZMs as a function of /i, when 
the normal state has a band inversion only at the T = 
(0, 0, 0)-point, i.e., m-p < 0. At [i — 0, the bulk is gapped 
and must have a pair of MZMs on opposite surfaces in a 
slab geometry. As /x is raised, these MZMs leak deeper 
into the bulk, but survive even after /i crosses |mr| de- 
spite the bulk now having a FS in the normal phase, 
gapped by superconductivity. A VPT eventually occurs 
at fi c = 0.9, at which the vortex is gapless and the sur- 
face MZMs merge into vortex line. Beyond /j, c , there are 
no longer any protected MZMs on the surface. 

Continuum limit: In the continuum limit of the lattice 
model, we can analytically calculate \i c . For k z = and 
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Figure 1. The vortex phase transition. Evolution of the low- 
est bulk/vortex states at k z = (top row), of the dispersion 
within the vortex (middle row) and of the surface MZMs (bot- 
tom row) as /i is varied when the normal phase has a band 
inversion at Y. At \i = 0, the normal phase is a strong TI 
and a superconducting vortex traps a MZM at its ends. As fi 
is increased, it first enters the conduction band at /i = |mr| 
and mid-gap states appear inside the vortex. For [i < fi c , the 
vortex stays gapped, but with a minigap S smaller than the 
bulk gap. The MZMs remain trapped near the surface. At 
/i c = 0.9, the gap vanishes signalling a phase transition. Be- 
yond jj, c , the vortex is gapped again, but there are no surface 
MZMs. We used the lattice Hamiltonian with the parameters 
t = 0.5, M = 2.5 and mo = —1.0. The pairing strength is 
Ao = 0.1 far away from the vortex and drops sharply to zero 
at the core. Other gap profiles give similar results. 



small k X:V around T, H k att reduces to the isotropic form 
T~Lk = vdt x <t • k + (m — ek 2 )r z — /i. In this form, a band 
inversion exists if me > 0. Thus, me < (> 0) defines a 
trivial insulator (strong TI). At k = v/m/e, mj. = m—ek 2 
vanishes and Hk resembles two copies of a TI surface. In 
particular, the Berry phase around each t x — ±1 FS is 
7T. We show later that this leads to a pair of vortex zero 
modes, signaling the VPT at fi c — vd\/ m/e. 

We solve analytically for the two bulk zero modes at /i 
to first order in A assuming |A(r)| = A 6(r — R), where 
Q is the step-function and ijlR/vd ^> 1. Calculating the 
zero modes separately for r < R and r > R and matching 
the solutions at r = R gives a pair of zero modes, only 
when fi = VD\Zm/e, for all vortex orientations. This is 
precisely where the momentum dependent 'mass' term 



changes sign [19]. Using the model parameters and the 
linearized approximation gives an estimate of /i c sa 1, in 
agreement with the lattice numerics. 

General Fermi surface Berry phase condition: 
For weak pairing, the VPT is expected to be governed 
by properties of the bulk FS. For concreteness, we begin 
by assuming we have a single FS in the k z = plane, 
which will be doubly degenerate due to the combined 
symmetry TI. We now argue that the VPT occurs when 
an appropriately defined Berry phase for each of the two 
degenerate bulk FSs is w. 

A convenient model for the vortex is A(r) = ^ — 
iy). The linear profile here simplifies calculations, but 
does not affect location of the zero mode. The choice of 
£ as the length scale gives the right minigap scale for the 
low energy excitations. Working in momentum space, we 
substitute r by id k , which gives 
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transforming now to the band basis | ip^) , which are eigen- 
states of the band Hamiltonian Hk\ipk) = E\cfk). Since 
we are only interested in very low energy phenomena, we 
project onto the two degenerate bands near the Fermi 
energy v = 1, 2. The projected Hamiltonian then is: 
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where D ka — d ka — iA a (k) and A a (k), the SU(2) con- 
nections, are 2x2 matrices: [A]^(fc) = i{Vk\dkJ\<p l C) ■ 

(i) Abelian case: Let us first consider the case when an 
additional quantum number (such as spin) can be used 
to label the degenerate FSs. Then, [A]^" must be diag- 
onal, and reduces to a pair of U(l) connections for the 
two FSs. In this situation, (3) is identical to the effective 
Hamiltonian for a p x + ip y SC, if we interpret momenta 
as position and ignore the gauge potential. The diagonal 
terms represent a transition from weak to strong pair- 
ing phase on crossing the FS when E k = [i [3]. Thus 
mid-gap are expected, composed of states near the Fermi 
energy. Due to the finite size of the FS, these states have 
an energy spacing of C(^|); the minigap energy scale. 
However, a zero energy state appears if the FS encloses 
a 7r-flux [3[. This can be implemented via the gauge po- 
tential if $ Fg A • dl = 7r leading to a pair of zero modes, 
since the other FS has the same Berry phase by T. 

(ii) General case, SU(2) connection: In the absence 
of any quantum number distinguishing the bands, one 
integrates the vector potential A(fc) around the FS in 
the k z = plane, to give the non-Abelian Berry phase: 
U B =Vexp [i§ ps A ■ dl] € SU(2), where V denotes path 
ordering. (There is no U(l) phase by T symmetry.) Al- 
though Ub itself depends on the choice of basis, its eigen- 
values e ±l ^ B are gauge invariant. A semiclassical anal- 
ysis [19] gives the Bohr-Sommerfield type quantization 
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condition for the low energy levels: 

E n = p-{2Tin + iT±^ B ) (4) 

where n is an integer and If is the FS perimeter. A pair 
of zero modes appears when (f>B = 7r, i.e. when Ub = 1 

We have considered a single closed FS in the k z = 
plane. Such a FS necessarily encloses a T invariant mo- 
mentum (TRIM), (e.g. T), given the symmetries. When 
there are multiple FSs, the condition above is applied in- 
dividually to each FS, since tunneling between them is 
neglected in the semiclassical approximation. Closed FSs 
that do not enclose a TRIM, or pairs of open FSs, cannot 
change the vortex topology. 

Candidate materials: We now apply the Berry 
phase criterion to some candidate materials to see which 
of them can have protected MZMs at the ends of vortices. 

Cu x Bi 2 Se 3 : The insulating phase of Bi 2 Se 3 is a strong 
TI with a single band inversion occurring at the T point. 
On Cu doping, Bi 2 Se 3 becomes n-type with an electron 
pocket at T and is reported to superconduct below T c — 
3.8 K [8, 20]. Photoemission measurements show fj, sa 
0.25 eV above the conduction band minimum at optimal 
doping [21]. We calculate the Berry phase eigenvalues 
for a FS around the T point numerically as a function 
of n, which evaluates to ±ir at pL c sa 0.24 eV above the 
conduction band minimum for a vortex along the c-axis 
of the crystal [19]. Hence /Lt > \i c indicates c-axis vortices 
are near the topological transition [19]. However, tilting 
the vortex away from the c-axis is found to raise /j, c to 
upto /i c = 0.30eK, for a vortex perpendicular to the c- 
axis. Therefore, sufficiently tilted vortices should host 
MZMs at the experimental doping level. 

p-doped TlBiTe 2 , p-doped Bi 2 Te 3 under pressure and 
Pd x Bi 2 Te 3 : The bands of TlBiTe 2 and Bi 2 Te 3 are topo- 
logically non-trivial because of a band inversion at the 
r point [22]. The topological character of Bi 2 Te 3 is be- 
lieved to be preserved under a pressure of up to 6.3 GPa, 
at which it undergoes a structural phase transition. On 
]>doping to a density of 6 x 10 20 cm~ 3 (3-6 x 10 18 cm -3 ), 
TlBiTe 2 (Bi 2 Te 3 under 3.1 GPa) becomes a SC below 
T c = 0.14 K (~ 3K) [10, 23], making it a natural sys- 
tem to search for the possibility of MZMs. Similarly, 
n-doping Bi 2 Te 3 to a concentration of 9 x 10 18 cm~ 3 by 
adding Pd reportedly results in T c = 5.5 K [8] in a small 
sample fraction. The superconductivity in Bi 2 Te 3 under 
pressure, and in TlBiTe 2 (Pda;Bi 2 Te 3 ) is believed to arise 
from six symmetry related hole (electron) pockets around 
the T-T line. This is an even number so vortex lines in 
superconducting TlBiTe 2 and both p- and n-type Bi 2 Tc 3 
should have MZMs at their ends in all orientations. 

MZMs from trivial insulators: The bulk criterion de- 
rived does not require a 'parent' topological band struc- 
ture. As a thought example, say we have four TRIMs 
with Hamiltonians like the continuum Hamiltonian "Hfe 
in their vicinity. Such band inversions at four TRIMs in 



a plane leads to a trivial insulator [24]. However, if their 
critical chemical potentials /i c differ, then there could be 
a range of fi where there are an odd number of VPTs be- 
low and above /i, leading to topologically non-trivial vor- 
tices. Interestingly PbTe and SnTe are both trivial insu- 
lators with band inversions relative to each other at the 
four equivalent L points. They both exhibit supercon- 
ductivity on doping below T c = 1.5 K [25] and 0.2 K [26] 
respectively. A combination of strain (to break the equiv- 
alence of the four L points) and doping could potentially 
create the scenario described above in one of these sys- 
tems. GeTe is similar to SnTe with T c ~ 0.3 K [27] but 
undergoes a spontaneous rhombohedral distortion result- 
ing in the desired symmetry. Thus, I- and T-symmetric 
systems with strong spin orbit can lead to SCs with 
vortex end MZMs, even in the absence of a proximate 
topological phase. Investigating the Fermi surface SU(2) 
Berry phases, and thus the vortex electronic structure, in 
this wide class of systems is a promising future direction 
in the hunt for Majorana fermions. 

In closing, we note that the VPT could potentially be 
probed via thermal transport along the vortex line. A 
hurdle is the small minigap scale (A/kp£, ~ A 2 /Ep), 
and the long confinement length of the MZMs to the 
surface, which may be ameliorated by considering strong 
coupling SCs or materials such as heavy fermions where 
E F is reduced. 

We thank A. M. Turner, J. H. Bardarson, A. Wray 
and C. L. Kane for insightful discussions, and NSF-DMR 
0645691 for funding. In parallel work, L. Fu, J. C. Y. Teo 
and C. L. Kane have arrived at similar conclusions. 
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CONTINUUM MODEL: ANALYTICAL 
SOLUTION 

In the continuum limit of the lattice model with cubic 
symmetry, we can analytically calculate \i c . For k z = 
and small k x , y around T, the lattice model reduces to 
the isotropic form ri k — vdt x <j • k + (m — ek 2 )r z — /i. 
In this form, a band inversion exists if me > 0. Thus, 
me > 0(< 0) characterizes a strong TI (trivial insulator). 
At k = y/m/e, mj. = to — ek 2 vanishes. Below we show 
that [i c = vpy/mfe is the critical chemical potential at 
which the VPT occurs. 

We solve for the two bulk zero modes at [i c analytically 
to first order in A by assuming |A(r)| = A o 0(r — R), 
where is the step-function R is a large radius such 
that R(A /vDfi) 3> 1. In the basis (tp(r), ia y ip^(r)) 
where ij){r) is a four-component spinor of electron an- 
nihilation operators indexed by spin and orbital indices, 
the Bogoliubov-de Gennes Hamiltonian in real space for 
a unit vortex is 
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where i/, are Pauli matrices in a combined orbital (rj) and 
spin ((Tj) space. For r < R, "H ra d has four zero modes 
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where R r = —ivr)T x cr • V + T z (m + eV 2 ) — fj,. The 0- 
dependence of "H BdG can be removed by observing that 
it commutes with the generalized angular momentum op- 
erator L z = —ids + a *t, n * ■ Since {L Z ,C} = 0, a single 
MZM, such as the one on the surface, must have an eigen- 
value n of L z , equal to zero. Now, the bulk zero modes 
can be thought of as the avenues through which the sur- 
face MZMs penetrate the bulk. Thus, they must have 
n = as well. The vortex preserves a mirror symmetry 
about the xy-plane described by M. = t z g z . Thus, R BdG 
can be block-diagonalized into sectors with opposite Ai 
eigenvalues. The two blocks are particle-hole conjugates 
since {M,C} = and each must contribute a single bulk 
MZM at fi c . The radial Hamiltonian for the M. = +1 
sector is 
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Matching the solutions and their derivatives at r = R to 
first order in Ao, using the asymptotic forms of the Bessel 
functions, is only possible at fj, = vr)\Jm/e. This gives 
an analytic solution to the critical chemcial potential [i c . 



where k± are roots of g(k) = e 2 k i +(l-2me)k 2 +(m 2 -(i 2 ) 
and Ji (x) are Bessel functions of the first kind of order 
i. For r > R, we drop all terms that contain r in the 
denominator. The zero modes of the remaining Hamil- 
tonian are of the form e tX±r x±, where A± is a root of 
e 2 A 4 + (1 - 2me)A 2 + (to 2 -(fi± iA ) 2 ) = with posi- 



VORTEX MODES AND THE BERRY PHASE OF 
THE FERMI SURFACE 

We consider a 3D system with s-wave pairing, where 
the bulk is pierced by a quantum flux h/2e. Assuming 
that the metallic Hamiltonian H k (the system without 



2 



superconductivity/quantum flux) has time-reversal sym- 
metry, we will derive the condition which governs the 
existence of zero vortex modes to the properties of the 
FS, namely the Berry phase. 

The mean field Bogoliubov-de Gennes (BdG) Hamil- 
tonian with s-wave pairing is of the form: 



(5) 



where = (c^, c^_ k (ia y )) is written in the Nambu ba- 
sis, capturing all orbital and spin degrees of freedom. The 
single particle Hamiltonian is 



H 



BdG, 



k,k') 



{H k - fi)5 kk > 
A*(r) (fi 



A(r) 



a y H -k a y) S kk' 



(6) 



where the pairing potential, given by A*(r) = 
(^I / |(r)^'|(r)), is position dependent due to the vortex. 

The quasi-particle operators ^ , are defined in such 
as a way that leads to spin singlet pairing. Since the 
metallic Hamiltonian H k has time-reversal symmetry, 
a y H*_ k a y = H k and the BdG Hamiltonian may be writ- 
ten as 



n 



BdG 



H k - y, A(r) 
A*(r) n-H k 



(7) 



Pairing potential 

In this section, we compute the pairing potential of a 
vortex in fc-space via a Fourier transform. We confine our 
system to be on a disk with radius £, the superconducting 
correlation length. The matrix element A k > k is 

A fe , fe = (fe'|A(r)|fe) 



1 

1 



d 2 re l(k - k,) - r A(r) 



rdrA{r) / dS e ^c OS (e-e q ) & - 

n JO 



where q = k — k' = q(cos6 q x + sin 6 q y). The 9 inte- 
gral evaluates to a Bessel function (of the first kind): 
2irie~ q Ji(qr). The matrix element becomes 



Afc'fc — 



2nie 



rdr A(r)J\(qr). 



(8) 



• At large q (q£ 3> 1), Ji(qr) sin qr/y/qr, and 
A(r) — > Ao becomes constant. The integral scales 
as ^-(qS,) 1 ^ 2 cosq£, and the matrix element A k ' k s» 
iA e- ie " (q£)~ 3/2 cos 9 £ ->■ 0. 

• At small g (q£ < 1), Ji(qr) ps gr/2 and the integral 
evaluates to (/Ao£ 3 /6. The matrix element scales 
as |Aoe~ 4e «q£ and is linear in g. 



(We drop the fc, fe' dependence, treating "H BdG as an op- 
erator.) 

We take the vortex to lie in the z direction, hence the 
pairing term takes the form A*(r) = A(r±)e , indepen- 



dent of z, where rj_e 



iO 



iy. The pairing ampli- 



tude becomes constant for large r±: A(r±_) — » Ao. Al- 
though the vortex breaks translational symmetry in the 
xy-plane, it is preserved in the z direction and hence, k z 
remains a good quantum number. The 3D system thus 
decouples into many 2D Hamiltonians enumerated by k z . 
Henceforth, we refer to r± as simply r, and k as the 2D 
momentum coordinate (k x ,k y ). 

The time-reversal operator is —ia y )C, taking k z — > 
—k Zl k — > —k, which is broken by the imaginary part 
of A. In addition, the system must have particle-hole 
(charge conjugation) symmetry given by the operator 
C = Tl y a y )C, where are the Pauli matrices acting on 
particle-hole space. Note that C also takes k z — > — k z , 
k — > —k. hence particle-hole is a symmetry of the 2D 
system only when k z ~ —k z (i.e. at or tt). 

In the remainder of this section, we do not assume 
anything about the value of k z , nor the symmetries of 
the 2D Hamiltonian H k \k z at a fixed k z . Our results 
remain valid applied to any 2D slice (with a smooth FS) 
of the 3D Brillouin zone, as long as the 3D Hamiltonian 
has time-reversal symmetry. 




Figure 1. The pairing potential. A fc / fc as a function of q£ 
where q = \k — fc'| and £ is the system radius. The wavefunc- 
tion tp(k) is smooth on the scale of £ , and hence the pairing 
potential may be modeled as —S'(q). 

We can model the pairing matrix element as a deriva- 
tive of the delta function: A k > k oc ie~ lBq {—5'{q)), this ap- 
proximation is valid in the regime where the wavefunction 
ijj(k) is smooth on the length scale of In Cartesian 
coordinates, the pairing term becomes iA e {dk x — idk y ), 
where A e is the effective p + ip pairing strength, with 
units energy x length -1 . From a simple analysis, we 
expect that A e sa Ao/£. 

The Hamiltonian in fc-space is of the form: 



H 



BdG 



H k 

*A e (<9 fcx 



— idh 



(9) 
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The eigenstates are of the form 

|V> = £w(fcM>, 



(10) 



where \<p£) is an eigenstate of Hk and v labels the band 
index. 



Projecting to the low energy states 

Solutions to the Hamiltonian (9) will consists of mostly 
states near the FS, where E k ~ /i. Hence we can simplify 
the system by projecting to the band eigenstates \ip%.) 
with energy near the chemical potential, where [i, v is 
the band index: 



an effective magnetic field V x A. Unlike a superconduc- 
tor, the effective field does not have to be localized nor 
quantized within the FS. 

Consequently, the spectrum of vortex states in our 2D 
system does not have any symmetry, and the zero modes 
are not topologically protected. We can only show that 
these modes are stable within the weak pairing limit un- 
der perturbation theory. We can restore particle-hole 
symmetry by combining the 2D systems at k z with that 
at — k z , at the cost of doubling the number of zero energy 
states. 

In the remainder of the section, we will explicitly show 
the following: when (an eigenvalue of) the Berry phase 
of the FS 4>f is n, there is an effective half quantum flux 
^ in the system which supports a Majorana mode [1-3]. 



(ii) 



If inversion symmetry is present, each band will be dou- 
bly degenerate. 

The projection of the metallic Hamiltonian Hk gives a 
diagonal matrix E£" = (<Pk\Hk\<p%) with its entries being 
the energies. The projection of the derivative operator 
gives the Berry connection: (<Pfc|«Vfc|</?£) = iVfc + A**", 
where the non-Abelian Berry connection is defined as: 
A"" = «K|V fe |^). Explicitly: 



H = 



A,. 



A ( i(d Px - id Py ) 
e \+A£ v - iA% v 



■ (12) 



The connections A x , A y are hermitian matrices which de- 
pends on the choice of Bloch states \<fk)- 

Dropping the band indices and writing Ef, as E(k), 
the Hamiltonian is: 



E(k) - /i 



A,. 



«(<9p* 
+A X 



+A X 



id Py ) 
- iA. 



id Py. 

■ iA„ 



-E(k) + n 



, (13) 



If we intepret momenta (fe) as position, then we can 
treat A x>y as a gauge potential and the Berry curvature 
F = dk x A y — dk y A x — i{A x ,A y ] as an effective magnetic 
field in the model. This system was studied by Read, 
Green, Ludwig, Bocquet and Zirnbauer in context of a 
Dirac Hamiltonian with random mass as well as p + ip 
superconducting systems [1-3]. 

While superficially similar to a p + ip superconductor, 
there is an important distinction - our system may not 
have particle-hole symmetry (unless k z = or w), and 
hence does not belong to the D class. To illustrate how 
C is broken, suppose that there is only one FS and hence 
Ai are simply lxl (Abelian) matrices. The system re- 
sembles a superconductor with p + ip pairing, but with 



Explicit solution with rotational symmetry 

This section is not necessary to the solution, but is in- 
structive and aids in the understanding of what the terms 
in the more general solutions mean. For simplicity, we 
only consider a single FS at wavevector fcp. We also as- 
sume an Abelian Berry connection, so A x ,A y are simply 
real numbers. 

With rotational symmetry, we can simplify the expres- 
sions in polar coordinates: k x +ik y = ke lB : E(k) = E(k), 
d kx -idk v =e- M (d k -f:de), A x -iA y = e^ b \A k - |A e ). 
In addition, it is possible to find a gauge for which 
Ak = 0, and Ag is a function of k, but independent of 9. 
Explicitly, the Berry connection is 

2nA (k) = 2Tr k'dk' F(k'), (14) 
Jo 

where F(k) = x A is the Berry curvature. The left 
side of (14) is the Berry phase along a circle of radius k, 
while the right side is the integrated Berry curvature. 
Our Hamiltonian simplifies to 



n 



iA e e 



E(k)- 
9 (d k + 



iA e e 



-i6 fa _ ide+Ag 
\ u k k 

-E(k) + /i 



(15) 



The Hamiltonian commutes with J z = —ide + H z , that 
is to say that the solutions are of the form 



J_ /w(fc)e i ("- 1 ) eN 
y/k V -iv{k)e in9 ' ' 



(16) 



for integers n (required by the wavefunction being single- 
valued) . Via the transformation 



-i(n-l)6 



i e 



-inO 



(17) 
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the effective Hamiltonian for (u, v) is 



n n = w ph nw- h L 



E(k) - fi 



i-h-At 



(18) 



Notice that the Hamiltonian is symmetric except for the 
terms proportional to dk, due to our choice of W p h- 

Our assumption is to replace k by the Fermi wavevector 
kp, Aq by Ag = Ag(kp) and ^ by This is justified as 
the amplitude of the wavefunction |u(fc)|, |v(fc)| is largest 
at the FS k — kp and exponentially decays away from 
the FS. The resulting Hamiltonian is equivalent to the 
Jackiw-Rebbi model [4]: 

U n = (n - A% - |) n x + iA e d k Tl y 

fx, p 1 

+ (E(k) - /i)n„. (19) 

A midgap state exists whenever E(k) — fx changes sign, 
with energy 



2tt 



(20) 



where <fip = 2irAg is the Berry phase of the FS. Hence, 
a zero energy solution exists when <pp is an odd multiple 

Of 7T. 

Explicitly, the eigenstates are of the form [1, 4] 



u(k) = exp / /3{k') dk', 



(21) 



where f3(k) must a be a decreasing function of k for u{k) 
to be normalizable. Assuming that E(k) is an increasing 
function of k, then 



(j, - E(k) 



v(k) = u(k) 



(22a) 
(22b) 



By inspection, this satisfies the Schrodinger equation for 
the Hamiltonian (19): 



Hr. 



f 



-A e £ n + A e /3' 
- A e /3 A e /3 



(23) 



For a Fermi velocity vp, E(k) — fx rj hvp(k — fc^). We 
can see that u(k) is a Gaussian with width WA~^]fVVp ~ 
^ Ao/hvpS,, which justifies the substitution k ^ kp ear- 
lier. 

We would like to point the reader to one last detail 
before moving on to the general solution. Putting our 
solution back in to the wavefunction (16) gives 



u{k)e v 



-i0 



(24) 




Figure 2. (Color online) Parameterization of the reciprocal 
space with orthogonal coordinates (E,rj). The (blue) closed 
contours are plotted for constant energies E, with the bold 
contour being the FS at the chemical potential /x. The (pur- 
ple) open contours are at constant 77, running from to 2n. At 
each point, t(E,r/) = g= and n(E,r]) — are vectors tan- 
gent and normal, respectively, to the energy contours. Notice 
that, for a closed FS, the direction of the vectors t, n rotate 
once counterclockwise as r\ is varied around the surface. 



The pseudo-spinor (in the Nambu basis) is an eigenstate 
of t ■ II, where t is a vector tangent to the FS at ik,p, &)■ 
The pseudo-spin locking to the the momentum gives rise 
to the 7r phase around the FS in the Hamiltonian (18) 
and Eq. (20). 



General solution without rotational symmetry 

The solution is similar in spirit to the case with circu- 
lar symmetry. Now, we label the momentum by energy 
contours (E, rf) instead of (r, 9) and demand constant E 
contours to be orthogonal to constant 77 contours. Similar 
to 9, rj is periodic with periodicity of 2tt. 

The idea of the derivation is as follows. We rewrite 
the Hamiltonian in a Jackiw-Rebbi form as a function E, 
perpendicular to the FS. Let T 1 , T 2 , T 3 be matrices which 
anticommute with each other, and square to the identity 
matrix. Then the differential equation 



r 1 (iA e 9 B ) + zr 2 (^- Ai ) + r 3 /(7 7 ) 



(25) 



has a bound state when E — fj, changes sign, the bound 
state is an eigenvector of T 3 and has energy f(rj) [4]. The 
remaining rj degree of freedom governs the existence of 
zero energy states, by requiring the wavefunction ip{E, rj) 
to be single valued. 

Define the vectors tangent and normal to the energy 
contours t = n = |g. Clearly, t _L n. (Fig. 2). 
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The derivatives and connections in (E, 77) coordinates are 
related to those in the Cartesian coordinates: 



t x t y 



9fe, 



and we define the 2x2 matrix to be J 1 

= j ( d A , ( A A = j 



d k , 



where J is the Jacobian matrix: 



A E 

An 



J = 



d(k x ,k y ) 



jE JT) 

Jy Jy^ 



(27) 



(28) 



From J 1 J = 1, we can see that ■ t = 0, 3 V ■ n = 0, 
hence J E || n, J n \\ t. 

First we rewrite the pairing term in terms of E and 77. 

*(3fc x - id ky ) 



A x - iA y = d E 



d v , where 



iA E ), 



d r , = i(J2-iJ])(d r ,-iA v ). 



(29a) 
(29b) 



It is always possible to find a gauge transformation which 
eliminates Ae near the FS. The transformation is of the 
form U\=V cxp [i J E A E (E')dE'] , where V is the path- 



(26) ordering operator. Since d E U\ = iA E lj\, 



U A {d E -iA E )U\=d E 



(30) 



This transformation will alter A v , since the derivative d v 
acts on In general, it impossible to make both A E 
and A n disappear. 

We make the substitution J E —iJ E = IJ^Ie - * 8 ", where 
8 n gives the direction of the normal vector n. As t is 
perpendicular to n, J£ - U]j = -i\3^\e~ ie ^, Eq. (29) 
becomes 



d E = ie- i9 »\J E \d E , 

d r) = e- ld "\r\{d 11 -iA r) ). 



(31a) 
(31b) 



At the moment, our Hamiltonian (13) is of the form: 



H = 



E - /j 

A e e ie » (i\3 E \d E - ndr, - iA v )) 



A e e- ie » (i^ldE + lJ^-iAr,)) 
-E + n 



(32) 



r 



The angle 9 n rotates by +2ir around a closed FS as 
77 varies from to 2ir. The e " phase in the off diago- 
nal terms of the Hamiltonian give rise to a tt phase in 
the eigenstates. We transform this phase away via the 
unitary transformation: 



such that 



where 



E-fi A e D- 
A e D+ n-E 



D~ = i\3 E \d E -i\3 r '\{id v + A v ) 

D + = i\3 E \d E + i\J ri \(idr, + A v ) + i| J 77 



dOn 

drj 



(33) 

(34) 

(35a) 
(35b) 



We have ignored the term dE@n from the assumption that 
the FS is smooth (no cusps). The term 9^6* n is extremely 
important as it will give us the n Berry phase shift. 

At any fixed value of 77, we write the Hamiltonian in 
the form of Eq. (25): 

U P hHU! h = (E- f i)r z + A e D x r x + A e D y T y , (36) 



where D x and D v are the symmetric and the antisym- 
metric parts of off the diagonal elements D ± . 



D a 



D y = 



D+ +D~ 
2 

D+ - D~ 
2i 



i\3 E \d E + i 



= \3 T >\ id., 



2 dr, ' 
2 drj 



A 



(37a) 
(37b) 



In the first expression (37a), there is an extra term 

^T^TST' wmcn we can absorb in to the energy deriva- 
tive via the transformation 



cxp 



/ 



9 n dE 



d E exp 



/ft, 



dE 



9e - 9r, (38) 



where g v (E) — J|je| l^T ■ This introduces a term d, n g, n 
in D y , but is irrelevant as g v is single- valued. (In the 
case with rotational symmetry, the term exp jg v is k 1 / 2 
in (17).) 

The (E dependent portion of the) solution to Jackiw- 
Rebbi Hamiltonian (36) is: 



ip(E, 77) ~ u(rj) exp 



J A 



E-n 

W\ 



dE 



(39) 



by treating D y as a constant (independent of E) on the 
energy scale of y/ A e \i E \. The effective Hamiltonian is 



U v = -A e D y . 



G 



Finally, we solve for D y {r})u{rf) = for a Majorana 
vortex bound state, subject to the constraint that u{rj) 
is single-valued. We can solve for u(rj) explicitly: 



u(rj) = Vexp 



ifdrj' [A v {r ] l ) + 



2 drj 



u(0). (40) 



Since 9 n winds by 2tt around a closed FS, §\d n 9 n = n 
gives an overall phase of— 1. We have u(2n) = — Ubu(O), 
where the Berry phase of the FS is 



Ub =V exp 



p2ir 

A n drj 



(41) 



Hence a solution exists for every —1 eigenvalue of Ub- 



Vortex bound states 

We can solve for the spectrum Caroli-de Gennes- 
Matricon bound states [5]. For an arbitrary energy 
£, the 'angular' portion of the wavefunction satisfies 
—A e D v u(r]) = £u(rj), equivalent to: 



-id v u(ri) = 



AWi +A " + m) E ,„ u "' > - (42) 



Note that | | ~ 1 = |t|. The solution for u(rj) is a the 
path-ordered exponential 



u(t]) = Vexp 

The full solution to the Hamiltonian (32) is: 
ip(E,ri) oc u(r)) exp 



2 dr] 



«(0). 



(43) 



—i 



(44) 



While A v is a hermitian matrix, £|t|/A e and d v 9 n are 
simply numbers, hence the integral may be evaluated sep- 
arately for each term. The integral §^ v \t\dr) is simply the 

perimeter of the FS lp. The integral M^ijMl is 
for a closed FS, and for an open FS (modulo 27r). The 
single-valued requirement u(27r) = u(0) (for a closed FS) 
becomes: 

u(0) = - exp (i^~£ ) U B u(0). (45) 

From this, we can calculate the allowed energies for an 
arbitrary Berry phase: 

E n - ^{2nn + TT-4>), (46) 
If 

for integers n, where e 1 ^ are the eigenvalues of Ub- Note 
that the eigenvalues always come in pairs ±0 due to the 
particle-hole symmetry of the system (considering both 
k z and — k z slices of the BZ). 

Numerical calculation of the Berry phase for Bi2Se3 

The insulating phase of Bi 2 Se 3 is a strong TI, with a 
band inversion occurring at the T point in its rhombohc- 
dral Brillouin zone. Cu x Bi2Se3 is reported to supercon- 
duct below T c — 3.8 K [6, 7]. Before the superconduct- 



ing transition, the carrier (electron) density is approxi- 
mately 2 x 10 20 cm~ 3 from Hall measurements [7]. Using 
the effective eight-band model from Ref . [8] , we estimate 
fi theoretically in this material from the carrier density 
to be w 0.4 eV relative to the conduction band bottom. 
However, photoemission measurements show fi w 0.25 eV 
above the bottom of the conduction band at the optimal 
doping (a; = 0.12) [9]. 

The conduction band minimum is at the T-point. For 
small carrier densities, we expect the FS to be centered 
around T. Using the same model from Ref. [8], we deter- 
mine fj, c for this material by calculating the Berry phase 
eigenvalues for a FS around the T point numerically as a 
function of fi (See Fig. 3). 

This calculation is done by discretizing the FS contour 
and making the pair of Bloch functions continuous and 
the Berry connection vanish at all but one points on this 
contour. More precisely, we parameterize the FS via rj £ 
[0, 2w] and compute the eigenstates \ f v {r])) for v = 1,2. 
The phase is chosen such that A% v — (ri)\d v \Lp v (r])) = 
for < rj < 2tt along the FS. 

In general, the Bloch functions will not be single- 
valued, i.e. |< y 3 !/ (27r)) ^ ^(O)). The unitary trans- 
formation required to rectify the discontinuity at this 
one point is precisely the non-Abelian phase: \ip^(0)} = 
\U bY v \<£ v {2ir)) . Because of time-reversal symmetry in 
the normal phase of Bi 2 Se 3 , the Berry phase Ub € SU{2) 
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has eigenvalues e ±% ^ which come in complex conjugate 
pairs. Fig. 3 shows the variation of the Berry phase eigen- 
values 6 function of /i for the FS of Bi2Se3 in the 
k c = plane, where the c-axis is normal to the layers. 
Here, /i is measured relative to the bottom of the con- 
duction band at the T point, which is where the band 
inversion occurs. Hence \x c 0.24 eV, which happens to 
be close to the chemical potential for these doped mate- 
rials according to ARPES. 
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Figure 3. Berry phase eigenvalues <f>B for Bi2Se3 as a func- 
tion of fi relative to the bottom of the conduction band. The 
eigenvalues appear in ± pairs and are only defined modulo 
2n. The phases <f>B are zero (modulo 2n) at both the conduc- 
tion band minimum and at energies far above the conduction 
band. Clearly, the cf>B = 7T at [i c — 0.24 eV, signalling a 
VPT. Photoemission measurements show /i w 0.25 eV ~ fi c 
above the conduction band minimum at optimal doping [9], 
and hence Cu doped Bi2Se3 is predicted to be near the vortex 
phase transition. 



